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ABSTRACT 

Using a geometry consisting of a hot central Comptonizing plasma surrounded by a thin accretion disk, we 
model the optical through hard X-ray spectral energy distributions of the type 1 Seyfert galaxies NGC 3516 and 
NGC 7469. As in the model proposed by Poutanen, Krolik, & Ryde for the X-ray binary Cygnus X-l and later 
applied to Seyfert galaxies by Zdziarski, Lubinski, & Smith, feedback between the radiation reprocessed by the 
disk and the thermal Comptonization emission from the hot central plasma plays a pivotal role in determining the 
X-ray spectrum, and as we show, the optical and ultraviolet spectra as well. Seemingly uncorrelated optical/UV 
and X-ray light curves, similar to those which have been observed from these objects can, in principle, be ex- 
plained by variations in the size, shape, and temperature of the Comptonizing plasma. Furthermore, by positing a 
disk mass accretion rate which satisfies a condition for global energy balance between the thermal Comptonization 
luminosity and the power available from accretion, one can predict the spectral properties of the heretofore poorly 
measured hard X-ray continuum above ** 50 keV in type 1 Seyfert galaxies. Conversely, forthcoming measure- 
ments of the hard X-ray continuum by more sensitive hard X-ray and soft 7 -ray telescopes, such as those aboard 
the International Gamma-Ray Astrophysics Laboratory ( INTEGRAL ) in conjunction with simultaneous optical, 
UV, and soft X-ray monitoring, will allow the mass accretion rates to be directly constrained for these sources in 
the context of this model. 

Subject headings: galaxies: active — galaxies: individual (NGC 3516, NGC 7469) — galaxies: Seyfert — 
X-rays: galaxies 


1. INTRODUCTION 

An important characteristic of type 1 Seyfert galaxies (here- 
after Seyfert Is) is the highly correlated variability of the con- 
tinuum flux across a wide range of optical and ultraviolet wave- 
lengths. The short relative lags between the various bands and 
the observed optical/UV color variations have been interpreted 
as evidence of thermal reprocessing by the accretion disk of 
higher energy emission from a single, relatively compact source 
(Krolik etal. 1991; Courvoisier & Clavel 1991; Collin-Souffrin 
1991). The principal candidate for this emission is the power- 
law-like X-ray continuum which is believed to be produced via 
thermal Comptonization in regions of hot plasma near the cen- 
tral black hole. Observations of relativistically broadened iron 
Ka emission lines and the so-called Compton reflection hump 
provide further evidence that a substantial amount of X-ray flux 
is intercepted and reprocessed by a thin accretion disk in these 
systems. 

Recent multiwaveband monitoring observations of Seyfert Is 
have shown that the optical/UV emission is not related to the 
observed X-ray flux in any simple way. In June-July 1996, a 
~ 30 day UV and X-ray monitoring campaign of NGC 7469 
by the International Ultraviolet Explorer (WE) and the Rossi 
X-ray Timing Experiment ( RXTE) showed that although the UV 
and 2-10 keV X-ray light curves displayed some rough simi- 
larities in variability amplitudes and time scales, they were not 
related in a fashion that would naively be expected from thermal 
reprocessing (Nandra et al. 1998 ). In particular, the UV max- 
ima preceded similar features in the X-rays by ~ 4 days, while 
the minima in both light curves occurred nearly simultaneously. 
A more recent analysis of these data revealed, however, that the 
UV light curve is well correlated with the X-ray spectral index 
as well as with the extrapolated X-ray flux in the 0 . 1-100 keV 


band. This suggests that there is a causal connection between 
the X-ray and UV emission in this object (Nandra et al. 2000). 
In April 1998, an intensive 3-day optical/UV/X-ray monitor- 
ing observation of NGC 3516 by the Hubble Space Telescope 
(HST) and RXTE showed optical fluxes which changed only 
by ~ 3% (peak-to-peak) while the 2-10 keV X-ray flux var- 
ied by ~ 60%. Furthermore, the shapes of the variations in 
the optical and X-ray light curves were only marginally sim- 
ilar. Nonetheless, the variability at optical wavelengths, from 
3590 A to 55 10 A, was highly correlated with essentially zero 
lag between the different continuum bands, consistent with the 
thermal reprocessing model (Edelson et al. 2000). 

Chiang & Blaes (2001a; hereafter Paper I) have shown that 
the 1998 optical and X-ray observations of NGC 3516 can 
largely be explained by a thermal Comptonization/disk repro- 
cessing model similar to that proposed by Poutanen, Krolik, & 
Ryde (1997) to account for the differences between the hard 
and soft spectral states in the Galactic X-ray binary Cygnus 
X-l. This model was also applied in a simplified form by 
Zdziarski, Lubinski, & Smith (1999; hereafter ZLS) to both X- 
ray binaries and Seyfert galaxies in order to explain an apparent 
correlation between the X-ray spectral index and the relative 
magnitude of the Compton reflection component. In Paper I, 
we found that inversely-related changes in the size and X-ray 
luminosity of the Comptonizing plasma can produce the rela- 
tively large variations in the 2-10 keV band while causing only 
much smaller variations in the disk thermal reprocessed emis- 
sion. With the addition of disk emission due to internal vis- 
cous dissipation, we found that the optical and X-ray variations 
could be roughly matched while simultaneously satisfying an 
energy balance condition relating the bolometric luminosity to 
the radiative power available through disk accretion. 
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In this paper, we refine the calculations for NGC 3516 pre- 
sented in Paper I by explicitly fitting the optical and X-ray con- 
tinuum data and producing broad band spectral energy distribu- 
tions (SEDs) for each observing epoch. In addition, we apply 
these calculations to the 1996 UV/X-ray monitoring observa- 
tions of NGC 7469 in order to provide a more complete phys- 
ical picture for the correlations found by Nandra et al. (2000). 
We proceed by discussing the relevant optical, UV, and X-ray 
observations in §2. In §3, we describe the thermal Comptoniza- 
tion and disk reprocessing model including our modifications of 
the geometry of the Comptonizing region and empirical color 
corrections for the disk emission. The results for NGC 3516 
and NGC 7469 are presented in §4; and we discuss these re- 
sults and conclude in §5. 

2. OBSERVATIONS 

For the 2-10 keV fluxes and spectral indices of NGC 3516, 
we have fit the archival RXTE/Proportional Counter Array 
(PCA) data for the April 1998 observations. Our X-ray spec- 
tral model consisted of an exponentially cut-off power-law plus 
a Compton reflection component (the pexrav model in XSPEC 
v 10.0), a Gaussian line for the Fe Ka emission line and neutral 
absorption for the foreground column density. The cut-off en- 
ergy for the exponential roll-over was set at E c = 120 keV which 
is sufficiently high so that it does not affect the inferred power- 
law index in the 2-10 keV band. We find somewhat softer spec- 
tral indices than Nandra et al. (1999) who analyzed the simul- 
taneous observations by the Advanced Satellite for Cosmology 
and Astrophysics ( ASCA ). This discrepancy is due to the fact 
that these authors did not include a Compton reflection com- 
ponent in their analysis (K. Nandra 2001, private communica- 
tion). The optical continuum variations have been taken from 
Fig. 2 of Edelson et al. (2000), and the mean optical fluxes 
at 55 10 A, 4235 A, and 3590 A are from their Table l. 2 Be- 
cause of the small aperture of the HST/Space Telescope Imag- 
ing Spectrograph, we do not apply a starlight correction. We do 
make an extinction correction using E(B-V) = 0.042 (Schlegel, 
Finkbeiner, & Davis 1998) and the reddening law of Cardelli et 
al. (1989). The flux measurements at 1360 A made by Edelson 
et al. were not simultaneous with the optical data, and as these 
authors discuss, the UV measurements are affected by detec- 
tor gain changes induced by thermal variations. Therefore, we 
restrict our analysis of the disk thermal emission to the optical 
bands. Since spectral fitting of the X-ray data entailed inte- 
gration times of ~ 30-70 ks (see Table Al), the optical fluxes 
used in the SED fits were averaged over the duration of each of 
these epochs. If the optical emission is due to disk reprocess- 
ing, it is expected to be produced at relatively large distances 
(r > 10 14 cm) from the central black hole. In accord with this, 
we have shifted the optical light curves by the measured peak 
lag of -0.21 days relative to the 2-10 X-ray light curve (Edel- 
son et al. 2000). The times of each observing epoch, the X-ray 
fit parameters, and the corrected and time-shifted optical fluxes 
are given in Table Al. 

The optical and UV continuum fluxes for NGC 
7469 were obtained from the AGN Watch web page 
(http : / /astronomy . ohio-state . edu/— agnwatch/) 
We use the 4865 A and 6962 A continuum light curves mea- 
sured by Collier et al. (1998) and the modified 1315 A light 
curves of Kriss et al. (2000). Welsh et al. (1998) compared 


the ground-based optical spectra of Collier et al. to contempo- 
raneous HST spectra and found that the ground-based spectra 
are substantially redder due to contamination from starburst 
activity in the host galaxy. In order to determine the continuum 
emission due to the Seyfert nucleus, we used the difference 
between the HST and the mean ground-based spectra to es- 
timate the starlight contribution. The simultaneous UV data 
were obtained using IUE (Wanders et al. 1997; Kriss et al. 
2000), and Welsh et al. (1998) showed that the IUE and HST 
fluxes at 1315 A were consistent over a half-day period during 
the 30-day monitoring campaign. Thus, despite the relatively 
large aperture used for the IUE spectra, we assume that the 
starlight contribution to be negligible at this wavelength. A 
reddening correction using E(B-V) = 0.12 has been applied to 
the optical and UV fluxes. This correction includes a Galac- 
tic component as well as absorption intrinsic to NGC 7469 
(Kriss et al. 2000). The mean corrected optical and UV fluxes 
ar e F1315 = 11.4, F4865 = 1.59, and F 6962 = 0.75, all in units of 
10“ 14 ergcm~ 2 s" 1 A"* 1 . The X-ray fluxes and spectral indices 
are taken from Nandra et al. (2000) who fit the RXTE/PCA 
data with the same model we used for NGC 3516. The optical 
and UV fluxes have been linearly interpolated from the AGN 
Watch data at the observing times of the X-ray data. The fluxes 
at 4865 A and 6962 A were time-shifted by the measured lags 
relative to the 1315 A light curve: -1.0 days for the 4865 A 
data and -1.5 days for the 6962 A data (Collier et al. 1998). 

3. THE MODEL: THERMAL COMPTONIZATION GEOMETRY, 
ENERGY BALANCE, AND THE DISK SPECTRUM 

The geometry we consider is similar to that proposed by 
Shapiro, Lightman, & Eardley (1976) in their two- temperature 
disk model of Cygnus X- 1 : a central, geometrically thick region 
of hot plasma surrounded by a relatively cold, thin accretion 
disk. In principle, the hot central region could arise because of 
a thermal instability as originally envisioned by Shapiro et al. 
(1976), or it could consist of an advection dominated accretion 
flow (ADAF) (e.g., Narayan & Yi 1995) or any one of its more 
recent variations (Blandford & Begelman 1999; Narayan, Igu- 
menshchev, & Abramowicz 2000; Quataert & Gruzinov 2000). 
Following Poutanen et al. (1997), our goal is to infer the prop- 
erties of the disk and Comptonizing plasma from the radiation 
physics and thereby provide constraints for the dynamical mod- 
els of the inner accretion flow which may help distinguish be- 
tween the various candidates. 

A key feature of the present calculations is the thermostatic 
feedback effect between the X-ray emission from the Comp- 
tonizing plasma and the soft photons which are produced by 
thermal reprocessing of these X-rays by the accretion disk 
(Haardt & Maraschi 1991; Stem et al. 1995). ZLS noted that 
the degree of feedback can be adjusted by altering the geom- 
etry of the hot plasma, and they invoked this mechanism to 
explain the apparent correlation between the X-ray spectral in- 
dex, T, and the relative magnitude, F, of the Compton reflec- 
tion component in Seyfert galaxies and X-ray binaries. Fol- 
lowing Poutanen et al. 1997 and similar results Zdziarski et al. 
(1998) obtained by applying a similar analysis to the X-ray bi- 
nary GX339-4, ZLS proposed that the region of Comptonizing 
plasma takes the form of a uniform, optically thin sphere cen- 
tered on the black hole. The X-ray spectral index variations oc- 
cur because the disk penetrates the plasma sphere by different 


2 For the spectral fitting, optical and U V wavelengths have been de-redshifted to the source rest frame using z = 0.0088 and c = 0.0163 for NGC 3516 and NGC 7469. 
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amounts due to changes in the sphere radius, and/or changes 
in the disk inner truncation radius, r^. Greater overlap of the 
disk by the sphere results in a larger fraction of the X-ray emis- 
sion being intercepted by the disk and also in a larger amount of 
reprocessed flux. This produces an increased Compton reflec- 
tion component as well as more soft disk photons entering the 
hot plasma. The additional disk flux provides enhanced cool- 
ing of the plasma by increasing Compton losses; and the equi- 
librium temperature of the plasma decreases, thus producing a 
softer X-ray spectrum. Using estimates for the fraction of ther- 
mally reprocessed emission intercepted by the plasma sphere 
as a function of the ratio AWn/rs and an empirical relationship 
between the Compton amplification factor, A, and the spectral 
index of the thermal Comptonization continuum (Beloborodov 
1999), 

r = 2.33(4 (1) 


ZLS derived an R-T relation which is similar to the observed 
correlation found in the X-ray spectra of Seyfert Is. 

The radial distribution of the flux incident upon a razor-thin 
disk due to isotropic emission from an optically thin sphere is 
given by 

Finc(r) = !6^7pG")’ (2) 


where L x is the total thermal Comptonization luminosity. The 
dimensionless function h(r/r s ) is Eq. 1 of ZLS. It has the 
asymptotic behavior h{r/r s » 1) -> (ir/4)(r,/r) 3 and is approx- 
imately flat for rjr s ^ 1 with h( 0) = 7r. In Paper I, we calculated 
the disk thermal spectrum by assuming local blackbody emis- 
sion and estimated the local disk temperature as being due to 
a combination of reprocessed flux and emission due to internal 
viscous dissipation: 


7disk(0 ~ 


(i-a)F iDC (r) + Fvi X (r) 


<?B 


1/4 
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where a = 0. 15 is the assumed disk albedo and a b is the Stefan- 
Boltzmann constant. As in Paper I, here we model the flux due 
to viscous dissipation as 


Evisc(r) — 


3 GMM 
87T r 3 



(4) 


where M is the black hole mass, M is the mass accretion rate, 
and n is the “effective” disk inner radius which we hereafter 
take to be n = r g = GM/c 2 , appropriate for a maximal Ken- 
metric. 

In Paper I, we showed that the near constant optical contin- 
uum of NGC 35 16 in the presence of large variations in the X- 
ray flux could be explained in this model if the thermal Comp- 
tonization luminosity were inversely proportional to the radius 
of the plasma sphere, L x oc rj 1 . In addition, we found that the 
disk inner truncation radius had to be fairly constant on time 
scales longer than the 3-day observation in order to avoid ex- 
cessive color variations in the optical and UV bands. However, 
if one assumes that the contribution to the disk thermal emission 
from internal viscous dissipation is negligible, then the amount 
of X-ray luminosity required to produce the observed optical 
continuum levels result in 2-10keV fluxes which exceed the 
observed values by a factor ~ 4. The addition of viscous dissi- 
pation in the disk can reduce the needed X-ray luminosity sub- 
stantially; and as it turns out, the amount of viscous dissipation 
required to reproduce the mean X-ray spectrum also satisfies an 
energy balance condition: 


■ Lv + L, 


disk 


where L disk is the luminosity of the disk thermal emission due 
to viscous dissipation assuming that the inner disk is cut-off at 
the truncation radius and Uatx is the total luminosity from 
viscous dissipation for an untruncated thin disk extending all 
the way down to r/ : 

r OO 

Ldisk,untr = 47T / rdrF wisc (r). (6) 

/ r mini r f 

The difference Luntr-Ldisk can be interpreted as the accretion 
power that would otherwise emerge as local blackbody emis- 
sion from the inner regions of the disk that is instead being used 
to heat the Comptonizing plasma. 

If one neglects local viscous dissipation, the intercepted seed 
photon luminosity is given by L s = L X L S , where L s is a function 
only of the quantity r^a/G (ZLS). Since the Compton ampli- 
fication factor is the ratio of the thermal Comptonization lumi- 
nosity and the seed photon luminosity, A H L x /L s = 1/1*, it too 
is a function only of r^nAs* In particular, its value monotoni- 
cally decreases as A*„iinA* increases. For r^n/rs = 1 and using 
Eq. 1, the equilibrium photon index is T = 1.55. Since the fitted 
X-ray spectral indices for NGC 3516 are T ~ 1.65, this implies 
that the plasma sphere must overlap disk so that r s > r^n for 
this model, assuming zero viscous dissipation. 

The inclusion of viscous dissipation has two effects on the 
amplification factor. First, as we note above, it reduces the X- 
ray luminosity required to produce the observed optical emis- 
sion. Second, it acts as a constant offset to the ZLS value of L s . 
Both of these effects reduce the amplification factor for a given 
ratio rminAj and cause the associated X-ray spectral index to 
be larger (Eq. 1). For fits to the NGC 3516 data presented in 
Paper I, the viscous dissipation needed to fit the X-ray fluxes 
was large enough to require that the radius of the Comptoniz- 
ing sphere had to be smaller than the disk inner radius. In the 
context of the plasma sphere geometry, this implies a large gap 
between the hot plasma and the inner edge of the thin disk. In 
order to avoid such a gap, we have modified the shape of the 
central hot plasma for cases where smaller disk covering fac- 
tors (i.e., smaller values of spectral index) are required. For 
r s < r m in, we take the shape of the Comptonizing plasma to be 
an oblate ellipsoid with the variable r s denoting its semi-minor 
axis. For r s > n , the Comptonizing region is spherical as be- 
fore. Figure A1 illustrates these two cases. 

For r s > rmin, the radial distribution of thermal Comptonized 
flux incident upon the disk is given by Eq. 2. For r s < r^n, it is 

Finc(r) = / sinacosada / Idcj), (7) 

™ z rmnr s J 0 J..^ 

where a is the angle an incoming ray makes with the disk plane, 
cj) is its azimuthal angle, and l = /(r, a,0;r 5 ) is the path length 
through the Comptonizing region (see Fig. Al). Explicit for- 
mulae for /, a ma x» and 0™* are given in the Appendix. 

The seed photon luminosity for the Comptonization process 
is the disk flux which enters the region of hot plasma. For our 
geometry, this is given by 

/ OO 

g(r)I(r)rdr, (8) 

.Tun 

where /(r) = [( 1 -a)F\ nc + F v isc ]/tt is the disk surface brightness, 
r- rjr mn , and we have defined 

pOi max 

g(r) = 2 / <f> mtx (a,r)sinacosada. (9) 

Jo 
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In Fig. A2, we plot the ratio of the seed photon luminosity to 
X-ray luminosity, L s /L Xy as a function of rs/r^n for the case 
where the disk emission is due entirely to thermal reprocessing. 
For comparison, we plot the results for the sphere+disk model 
described by ZLS and the hybrid sphere/ellipsoid model. 

We include an empirical color correction for the disk ther- 
mal emission as a function of radius. Hubeny et al. (2000, 
2001) have modeled disk atmospheres appropriate for Seyfert 
nuclei and find that the optical continuum can be moderately 
well approximated by local blackbody emission, but at UV 
wavelengths and higher, significant deviations occur. This also 
seems to be true for disks powered solely by external illumina- 
tion (Sincell & Krolik 1997). Since we are modeling the UV 
emission for NGC 7469, a color correction, in principle, should 
be applied. We adopt a temperature dependent color correction 
of the form 


Tcol(Tdisk) — /oo 


(foe ~ 1)(1 +exp(-f/fr/Ai/)) 
l+GXptfj/p-I/frVAp) 


(10) 


where v p = l^lksT^/h is the peak frequency of a blackbody 
with temperature 7^. This expression for / co i goes from unity 
at low temperatures to / at high temperatures with a transition 
at v b » v p . We find that = 2.3 and v b = Av = 5 x 10 15 Hz 
do a reasonable job of reproducing the model disk spectra of 
Hubeny et al. (2001). 

We fit the observed optical, UV and X-ray data in the follow- 
ing manner: For a given trial value of plasma temperature, 7, 
the thermal Comptonization continuum is modeled as an expo- 
nentially cut-off power-law: 


L E = (Lo/£,)(£/£ f )'- r exp(-E/£ f ), (11) 


where we set E c = 2kgT appropriate for plasma Thomson depths 
r < 1. The normalization Lq is given by 


r rlOkeV/Ec / E \ HT 

, t y w* < i2 > 

where F 2-10 is the measured 2-10 keV flux and d\ is the lumi- 
nosity distance to the source. The total thermal Comptonization 
luminosity is then 


L x = Lo f exp (-E/E c )d(E/E c \ (13) 

7 E™ /E c \*W 


where we have set the minimum photon energy at a nominal 
value of = 0.01 keV. 

Having determined the thermal Comptonization luminosity 
for a given plasma temperature, we find the value of r s which 
provides the best-fit disk spectrum: 


4 ttcosi hv 3 f°° rdr 

d f c2 ft o, Jr mz exp(/ii/// col * fl 7' disk ) - r ( 1 4 ; 

where i is the observer inclination, and the disk temperature is 
given by Eq. 3. We then use the best-fit value of r s to compute 
the corresponding seed photon luminosity, L s , the Compton am- 
plification factor A - L x /L s , and the predicted value of the X-ray 
photon spectral index, r mode i. In contrast to ZLS and Paper I, 
we use the relation 


(15) 


rather than Eq. 1 of Beloborodov (1999). Malzac, Beloborodov, 
& Poutanen (2001) obtained this parameterization for pill-box 
shaped coronae atop an accretion disk using a Monte Carlo 
method to compute the thermal Comptonization continuum, 
and we have found a nearly identical relation from our own 


Monte Carlo simulations of thermal Comptonization and disk 
thermal reprocessing in the sphere+disk geometry (Chiang & 
Blaes 2001b). The best-fit value of the plasma temperature 7 
has been ascertained when the predicted photon spectral index 
matches the measured spectral index. The effective Thomson 
depth, T c fT, of the Comptonizing plasma is estimated by invert- 
ing another empirical relation found by Beloborodov (1999), 



where the Compton y-parameter is y = 49 e (\ +40 e )r c ff(l +r eff ) 
and 9 e = k B T Although the specific values of r eff that 
we find have no substantive role in our calculations, they may 
provide, along with the plasma temperatures, useful constraints 
on the properties of the Comptonizing region for the dynamical 
models. 

The other important model parameters are the disk inner 
truncation radius, r^n, the central black hole mass, Af, and the 
mass accretion rate, M . We assume that these quantities are 
constant throughout the monitoring campaigns for each object. 
The black hole masses can certainly be taken to be constant, 
and rmin and M will likely only change appreciably on viscous 
time scales which are on the order of O(K) 1 ) years for typical 
Seyfert 1 parameters. For NGC 7469, we have chosen a black 
hole mass of M 7 = (M/ 10 7 M©) = 1 . This is consistent with the 
upper range of the mass estimates by Collier et al. ( 1998) using 
broad emission line reverberation mapping. As we discussed 
in Paper I, the disk inner truncation radius will be largely con- 
strained by the shape and magnitude of the optical/UV contin- 
uum. For NGC 7469, the optical and UV fluxes show evidence 
for the roll-over associated with the peak of the black body 
spectrum at this inner radius. Since the temperature distribution 
of the disk will nearly follow a T&& oc r~ 3 ' 4 radial dependence 
whether the local disk flux is due to viscous dissipation or ther- 
mal reprocessing from an extended central X-ray source (see 
Paper I), we can use the optical and UV colors and magnitudes 
for NGC 7469 to estimate values for the disk inner radius 
and the disk temperature 7di S k,max at that radius. Insofar as the 
radial temperature distribution obtained from our more detailed 
thermal reprocessing plus viscous dissipation calculation does 
obey a r~ 3 ' 4 dependence, the values we find for /-„„„ and T^max 
will depend only on the optical and UV fluxes. Performing 
these fits over the 30 epochs of optical/UV data for NGC 7469, 
the average values we obtain are n = 3.0 ± 0.3 x 10 14 cm and 
7disk,max = 3.5 ± 0.2 x 10 4 K. As we showed in Paper I, disk 
spectra computed for values of r^n which differ significantly 
from the above will have very different optical/UV colors if 
normalized to the longest wavelength optical flux. As a poste- 
riori confirmation of setting r min = 3x 10 14 cm for NGC 7469, 
we note the quality of the fits to the optical/UV data shown in 
Fig. A7 and that the temperature distributions for those calcu- 
lations do indeed have close to a r -3 / 4 dependence. Now, since 
r nun £ 200 r^, the amount of disk emission due to viscous dis- 
sipation is insensitive to the term containing r/ in Eq. 4 and 
depends almost entirely on the product MM. Therefore, the full 
range of black hole masses, M ~ 10 6 — 1 0 7 M©, found by Collier 
et al. can be accommodated by adjusting M accordingly. 

Simultaneous ASG4 measurements of the broad Fe Ka line 
emission in NGC 3516 indicate that much of this emission orig- 
inates from disk radii < 6r ? , very close to the central black 
hole (Nandra et al. 1999). Therefore we adopt a disk inner ra- 
dius of rmm = 6 r g , so that the choice of black hole mass also 
sets the value of r^n- Unlike NGC 7469, the optical and UV 
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continuum fluxes we have for NGC 3516 do not show a high 
frequency roll-over, and we can infer only an upper limit for 
the disk truncation radius from these data. After de-reddening 
the optical and UV fluxes, a disk inner truncation radius of 
rmin < 1 x 10 13 cm seems to be required, yielding a black hole 
mass of M < 1O 7 M 0 . 3 However, a mass this small yields 
somewhat large values for the temperature of the Comptonizing 
plasma during some epochs (see Fig. A3), so we have computed 
models for several black hole masses Mi = 1, 2, and 3. 

Given values for n and Af , several considerations limit the 
plausible range of accretion rates. First, since they are assumed 
to be constant, they must be sufficiently small to allow the ob- 
served optical and UV variability to be due to disk thermal re- 
processing. Second, since the amount of disk flux and the ratio 
r s /r min determine the X-ray luminosity through the observed 
photon index and amplification factor scaling relation, larger 
accretion rates correspond to smaller thermal Comptonization 
roll-over energies, E c = 2 kgT. The values of E c cannot be too 
small, otherwise the roll-over would be evident in the observed 
X-ray spectra. Therefore, somewhat arbitrarily, we selected ac- 
cretion rates which ensured that kgT > 50 keV in all epochs. As 
it turns out, this condition is more stringent than the constraint 
that the varying part of the optical/UV flux need be due to ther- 
mal reprocessing. On the other hand, the accretion rate must 
be sufficiently large so that energy balance is satisfied such that 
Unti > £jc + £disk- Furthermore, larger accretion rates are pre- 
ferred since they imply values of kgT which are not too large. 
This is essential since we do not consider the effects of pair bal- 
ance or emission from non-thermal particle distributions, and 
thus this model tacitly assumes kgT < m e c 2 . The various pa- 
rameters which are held constant are summarized in Tables A2 
and A3. 

4. RESULTS 

4.1. NGC 3516 

The model parameters and resulting light curves for NGC 
3516 are shown in Fig. A3. For the Mi = 2 and Mi = 3 cases, the 
fitted plasma temperatures and optical depths are well within 
the range of values which have been typically measured or 
assumed for type 1 Seyfert galaxies (Zdziarski, Poutanen, & 
Johnson 2000). However, for the Mi = 1 case, the plasma tem- 
perature is somewhat high in two epochs, kgT ~ 1 keV. In all 
three cases, the size parameter of the Comptonizing region, r 5 , 
is about a factor of 3-4 smaller than the disk inner truncation 
radius so that the shape of the region containing the hot plasma 
is relatively flat. By comparison, the range of spherical radii 
we found for the fits presented in Paper I were only somewhat 
smaller than the truncation radius, r s ~ 0.7 r m j n . The reason for 
this less pronounced difference can be seen in Fig. A2 in which 
the L s fL x curve for the ZLS calculation declines more rapidly 
as r s /r m j n decreases for r 5 /r mn < 1 than it does for the hybrid 
sphere/ellipsoid model. In the lower right panel of Fig. A3, we 
see that the thermal Comptonization luminosities, L x , are con- 
sistently less than the available accretion power, Luna- - Ed, S k ( see 
Table A2). For the black hole masses and accretion rates we 
consider, the thermal Comptonization luminosity constitutes a 
substantial fraction of ntr -E d j s k, implying fairly efficient ac- 
cretion, comparable to that for an untruncated thin disk in a 
Kerr metric. 

The spectral energy distributions for these fits are shown 
in Fig. A4. The optical and X-ray spectra are reasonably 

In Paper I, the effects of reddening were not considered, and a truncation radius 


well described by the model SEDs, but the continuum flux at 
13 60 A is problematic. When the reddening correction is ap- 
plied, this datum sits above a vF v oc i/ 4/3 extrapolation from 
the optical data. Even if we made the model disk spectra bluer 
by, for example, decreasing the inner truncation radius below 
/“min = 8.9 x 10 12 cm, the 1360A flux would still exceed the 
model value by ~ 20%. Edelson et al. (2000) estimate sys- 
tematic changes in the UV flux of 1.5% due to thermal vari- 
ations, but this is too small to account for the potential dis- 
crepancy. Furthermore, given the lack of a strong X-ray flare 
either during or preceding the UV measurements (see Fig. 1, 
Edelson et al. 2000) and the ±1.5% relative variations in the 
optical light curves, which directly followed them, it is unlikely 
that there was a sufficiently strong flare in the disk emission 
to account for such a large UV continuum flux. However, the 
1360A continuum may be contaminated by the broad wings of 
the Si IV 1397A and Cll 1 335 A emission lines. It would be 
extremely useful for a more complete analysis of these data to 
be performed that is similar to the spectral fitting by Kriss et al. 
(2000) for the HST/Faint Object Spectrograph and 1UE spectra 
of NGC 7469 in order to determine the separate continuum and 
emission line fluxes as precisely as possible. 

4.2. NGC 7469 

The corresponding fits for NGC 7469 shown in Fig. A5 
(solid diamonds) present a different picture for the inner re- 
gions of this object and seem to require additional elements in 
the model. First of all, the generally softer X-ray spectra imply 
values of r s which are close to the disk inner radius so that the 
central plasma region is spherical, save for one epoch when the 
X-ray spectrum is its hardest. Of course, the specific values we 
fit for r s also depend on our choice of M. Second, the variations 
in the observed optical fluxes are noticeably smaller than that 
of the model. Part of the excess variability in the model opti- 
cal light curves is due to the fact that light travel time effects 
are not included in these calculations. The measured time lags 
of ~ 1 day imply that the model curves should be smeared out 
on at least this time scale. However, even if such light travel 
time effects are included, there will still be excess variability at 
the ~ 10 day time scale. Nonetheless, given the simplicity of 
this model, it is encouraging that the shape and magnitude of 
the UV light curve can be so well reproduced while the relative 
magnitude of the optical fluxes are approximately matched and 
the optical variability is not greatly over-predicted. 

Unfortunately, the relatively large variations in the UV flux 
seem to require much larger changes in the X-ray luminosity 
than are found by simply extrapolating the 2-10 keV flux over 
even as broad a range as 0. 1 to 100 keV (see Nandra et al. 2000). 
Given the measured 2-10 keV fluxes and photon indices, the 
size of the L x variations which are required by thermal repro- 
cessing can only be accommodated by setting the thermal roll- 
over energy to unrealistically large values — in several epochs, 
we find kgT > 10 5 keV. Clearly, such large temperatures are 
not consistent with the underlying assumptions of these calcu- 
lations. Not only would pair-production and non-thermal emis- 
sion likely be present at these temperatures (Coppi 1999), the 
shape of the thermal Comptonization continuum at 2-10 keV 
would not resemble a power-law as individual scattering or- 
ders would become apparent in these spectra for temperatures 
> 10 3 keV (e.g.. Stern et al. 1995). 

r iUin ^3 x 10 13 cm was sufficient 
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One way of narrowing the required temperature range is by 
giving the observed X-ray continuum more “leverage’* in pro- 
ducing a greater amount of variability in the thermal repro- 
cessed disk emission. This can be accomplished by allowing 
for anisotropy in the thermal Compton emission so that the X- 
ray intensity that we see directly is less than the X-ray intensity 
incident upon the accretion disk. Such anisotropy may arise 
in several different ways. The strong gravitational field of the 
central black hole will cause radiation produced near it to fol- 
low curved geodesics which may intercept the disk rather than 
follow straight-line trajectories to infinity (e.g., Cunningham 
1975). Bulk relativistic motion of the Comptonizing plasma 
will beam the radiation along the direction of motion (Be- 
loborodov 1999). However, in this case, the bulk motion would 
have to be directed towards the disk instead of away from it as 
would occur in a mildly relativistic jet. Finally, since the region 
of Comptonizing plasma is not being illuminated uniformly, but 
rather by the seed photons from the disk encircling its equator, 
there is a strong preference for the thermal Compton photons of 
the first scattering order to be directed back towards the disk. 

Ideally, we would have full descriptions of the geodesics fol- 
lowed by the thermal Compton radiation in order to account 
for the aforementioned effects properly, but the use of a simple 
parameterization will serve to illustrate how large a difference 
some mild anisotropy can make. Therefore, rather than model- 
ing the various effects in detail, we instead introduce a param- 
eter, which describes the anisotropy so that the total lumi- 
nosity of the Comptonizing plasma and the luminosity used to 
compute the thermal reprocessed emission is given by 

= (17) 

Here L Xy app is the apparent luminosity inferred from the assumed 
thermal roll-over, E c = 2 k B T, the measured spectral index, and 
the measured 2-10 keV flux, i.e., it is the L x given by Eqs. 12 
and 13. Setting £ = 1.5, we have recomputed the spectral fits 
to the NGC 7469 data and plotted the results in Fig. A5 (plus- 
signs). As expected, a narrower range of plasma temperatures 
are found which reproduce the optical/UV variability, but the 
plasma temperatures nonetheless reach values k B T ~ 10 4 keV, 
which are still too large to be self-consistent. 

We now ease some of the restrictions on the parameters of the 
X-ray spectrum. For all previous fits, we have used the best-fit 
values of the photon indices obtained from the X-ray spectral 
fitting. Since many of the spectral indices for the NGC 7469 
observations are near T & 2 (see Fig. A6), a given change in 
the thermal roll-over energy E c has less impact on the overall 
X-ray luminosity than it would for a harder spectrum with, say, 
T « 1.6. Therefore, in order to increase the accessible range 
of the X-ray luminosities, we allow the X-ray spectral index to 
vary within the 1-sigma bounds of the best-fit values, and we 
allow for small variations in the 2-10 keV flux used to compute 
L x (Eq. 12). Furthermore, we restrict the plasma temperature to 
be 40 keV <k B T < 700 keV; and we set the anisotropy param- 
eter to be f = 1.5, as before. With these additional conditions, 
we obtain the spectral fits shown in Fig. A6. For some epochs, 
the X-ray spectral indices attain values which are right at the 
1 -sigma limits, and small deviations from the observed X-ray 
fluxes and somewhat larger deviations in UV fluxes are now 
present, but the overall fits of the spectral energy distributions 
are largely satisfactory while maintaining a consistent range of 
plasma temperatures. We plot in Fig. A7, the SEDs for six ex- 
emplary epochs of the NGC 7469 data for the three cases we 


have considered. The details of these three fitting schemes are 
given in Table A3. 

5. DISCUSSION AND CONCLUSIONS 

Given the parameters from the above spectral fits, NGC 3516 
and NGC 7469 provide contrasting cases for the characteristics 
of the accretion flow in the inner regions of Seyfert galaxies. 
NGC 3516 has a high accretion efficiency (> 30%), a small 
disk inner truncation radius, and a flattened, disk-like distribu- 
tion of Comptonizing plasma. On the other hand, NGC 7469 
has a substantially lower accretion efficiency (~ 2-3%), a large 
truncation radius, and a vertically extended hot plasma region. 
As we have noted above, the small truncation radius for NGC 
3516 is consistent with the relativistically broadened iron K a 
line measured by ASCA over the same observing period. Al- 
though similar ASCA observations were not taken during the 
1996 campaign for NGC 7469, Nandra et al. (2000) found iron 
K q equivalent widths of ~ 150eV from the RXTE/PCA data. 
This value agrees with the equivalent width measured during 
the November 1993 ASCA observations by Guainazzi et al. 
(1994) who found the iron Ka emission line to be very nar- 
row. If the Fe Ka line is composed of emission from a disk, 
then Guainazzi et al. concluded that it must originate at least 
several tens of Schwarzschild radii from the central black hole, 
consistent with our estimate for the disk inner truncation radius 
of r mn ~ 200 r^. 

As we discussed in the Introduction, a similar analysis was 
performed by Poutanen et al. (1997) for the X-ray binary 
Cygnus X-l. With respect to the semi-analytic model they used, 
the main difference between that work and the present is in the 
empirical relation employed by those authors to characterize the 
thermal Comptonization spectrum in terms of the model param- 
eters. In particular, rather than Eq. 15, Poutanen et al. used the 
scaling relation 



which was found from Comptonization modeling by Pietrini & 
Krolik (1995). 4 On the observational side, since Cygnus X- 
1 is significantly brighter in the hard X-rays than any Seyfert 
galaxy, both the roll-over energy of the thermal Comptoniza- 
tion continuum and the amplitude of the Compton reflection 
component for Cygnus X-l could be measured directly. These 
data provided Poutanen et al. with estimates of the plasma tem- 
perature as well as the ratio r 5 /r^ n . In addition, the thermal 
disk spectrum, which was observable in the soft X-rays, gave 
them constraints on the accretion disk flux and the inner disk 
radius. Armed with this information, Poutanen et al. were able 
to apply the scaling relation and the assumed sphere+disk ge- 
ometry to determine the mass accretion rate directly in both the 
hard and soft spectral states rather than to have to assume values 
for the accretion rates based on “reasonable” ranges for param- 
eters such as the plasma temperature, as we have done here. 

As the analysis of the Cygnus X- 1 data demonstrates, more 
stringent tests of this model in its application to Seyfert galax- 
ies will be available when the hard X-ray/soft gamma- ray spec- 
tra from these objects can be reliably measured on sufficiently 
short time scales to detect the changes in the spectral roll-over 
implied by the temperature curves in Figs. A3, A5, and A6. 
This will require a hard X-ray telescope with a large effec- 
tive area at these energies such as that proposed for ASTRO- 


4 Eq 18 differs numerically from the expression used by ZLS, Eq. 1, by only < 3% for photon indices ranging from T - 1 .4 to 2.3. 
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E2. Nonetheless, relevant constraints for the mean properties 
of the Comptonizing plasma may be ascertained from exist- 
ing longer time scale observations by the Beppo-SAX satellite 
and the Oriented Scintillation Spectrometer Experiment (OSSE) 
aboard the Compton Observatory and from forthcoming ob- 
servations which will be available with INTEGRAL. For ex- 
ample, spectral fits of data from OSSE observations of NGC 
4151 (Johnson et al. 1997) and from Beppo-SAX observations 
of NGC 5548 (Nicastro et al. 2000) suggest that as the X-ray 
power-law spectrum softens, the cut-off energy (or equivalently, 
the plasma temperature) increases, contrary to naive expecta- 
tions from thermal Comptonization models. However, since the 
Thomson depth of the plasma may also change simultaneously, 
a sufficiently large decrease may more than compensate for any 
spectral hardening which a temperature increase would other- 
wise imply. In Fig. A8, we plot the plasma temperatures versus 
the measured X-ray spectral indices for the “optimal” case fits 
of the NGC 7469 data (i.e., £ = 1.5 and T allowed to vary) and 
for the three cases we considered for NGC 35 16. It is clear that 
a strict relationship between k B T and T is not required. Since 
changes in k B T will affect both L x and the thermal reprocess- 
ing contribution to the disk emission, simultaneous optical and 
UV monitoring will be necessary in order for longer term mea- 
surements of the hard X-ray spectra of Seyfert Is to constrain 
thermal Comptonization models in which thermal reprocessing 
plays an important role. 

Greater spectral coverage in the hard X-ray band will also 
provide improved constraints on the apparent correlation be- 
tween the thermal Comptonization spectral index and the rela- 
tive strength of the Compton reflection component. However, 
the predictions of the ZLS analysis are altered since the inclu- 
sion of disk emission due to viscous dissipation implies quite 
different estimates for the seed photon luminosities. The X- 
ray spectral indices will no longer be solely a function of the 
geometry, and simultaneous optical and UV monitoring along 
with broad band X-ray spectroscopy will again be necessary 
to determine the strength of the Compton reflection component 
implied by this model. 

In fact, for the present data, we do find significant differ- 
ences between the predicted reflection fractions and the values 
obtained by fitting the PEXRAV model to the RXTE/PCA data. 
In addition, the reflection fraction given by this model for a 
specific viewing angle can differ substantially from the aver- 
age Compton reflection strength as defined by ZLS. Following 
ZLS, the total thermal Comptonization luminosity intercepted 
by the disk is 

r OO 

L d =4iT drrF inc (r), (19) 

f \ min 

where the incident flux Fj nc is given by Eqs. 2 or 7. This is the 
luminosity intercepted by both sides of the disk. ZLS define 
the reflection fraction, averaged over all viewing angles, to be 
the ratio of the luminosity intercepted by the disk and the ther- 
mal Comptonization luminosity which can be seen by distant 
observers: 

RzLS = r ^ r - ( 20 ) 

When comparing the observed R-T correlation for a large sam- 
ple of objects which have random orientations, this is the ap- 
propriate form to use for the reflection fraction. However, for 
a specific object with a known inclination, this expression may 
not be a good approximation. In particular, for an object with 
disk inclination i = 0 and plasma radius r s < r min , the full X-ray 


luminosity from the Comptonizing plasma is seen, but only half 
of the disk reprocessing emission is observed. In this case, one 
should use 

Rh o = y^l- (2D 

Lf.app 

Note that the apparent X-ray luminosity appears in the denomi- 
nator, but Ld in the numerator is computed using L x = £L Xiapp 
to account for any anisotropy. In Fig. A9, we plot the R- 
T relationships using these two expressions for the reflection 
fraction along with the fitted values from our analysis of the 
RXTE/PCA data for NGC 3516 and from the analysis of Nan- 
dra et al. (2000) for NGC 7469. The model values using R i=0 
clearly under-predict the fitted values by a substantial amount 
and are also much smaller than the average values, Rzls As 
Nandra et al. (2000) point out, any apparent correlation be- 
tween spectral index and Compton reflection strength may be 
partially due to a statistical correlation between the parameters 
which arise from model fitting rather than from intrinsic prop- 
erties of the source. The presence of such a statistical correla- 
tion between fitting parameters suggests that the small reflec- 
tion fractions given by this model may be accommodated by 
the data if the underlying thermal Comptonization spectra are 
actually somewhat harder than those found by fitting the simple 
reflection models to the RXTE/PCA data. Since RXTE obser- 
vations generally do not provide much useful spectral coverage 
above ~ 15 keV for Seyfert galaxies, more sensitive hard X-ray 
continuum measurements for energies >10 keV are required to 
resolve this discrepancy. 

Even in the absence of high quality hard X-ray data to mea- 
sure the thermal roll-over or the Compton reflection fraction 
definitively, certain aspects of this model are indeed tested by 
the available data, though the constraints are not as direct. For 
example, if the UV fluxes for NGC 7469 were consistent with 
an extrapolation of the vF v ~ dependence of the optical 
data, then a much smaller disk inner truncation radius would 
be required. This in turn would imply a much broader Fe Ka 
emission line that may contradict the observations of a nar- 
row line by Guainazzi et al. (1994), or it may require a much 
smaller black hole mass, which could be in conflict with the 
reverberation mapping estimates of Collier et al. (1998). Fur- 
thermore, the quantity MM ~ 10 6 M 0 2 yr 1 which we find for 
NGC 7469 (see Table A3) is in rough accord with the value of 
MM ~ 0.7 x 10 6 Mq 2 yr -1 found by Collier et al. (1998) by 
fitting interband continuum lags for these same observations. 

The large disk inner truncation radius and vertically extended 
Comptonizing region we find for NGC 7469 suggests that the 
inner regions of this object may consist of an advection domi- 
nated accretion flow (ADAF) similar to those which have been 
used to model the X-ray spectra of Galactic X-ray binaries (e.g., 
Esin, McClintock, & Narayan 1997; Esin et al. 2001). How- 
ever, the mass accretion rate which we find for NGC 7469 in 
our optimal case, M ~ O.15M 0 yr -1 , is larger than the Ed- 
dington limit by a factor of ~ 3, assuming a Kerr metric and 
a 10 7 M e black hole; and it is - 40% of the Eddington limit 
for a standard thin disk in a Schwarzschild metric. This ac- 
cretion rate likely exceeds the critical accretion rate found by 
Narayan & Yi (1995) of Af cri t ~ a 2 Af E dd, below which advection 
dominated flow is expected to occur. The anomalous viscos- 
ity parameter is found in magnetohydrodynamic simulations of 
disks to have values a < 0. 1 (e.g., Balbus & Hawley 1998). As 
we noted above, we chose the largest possible accretion rates 
for our models which were consistent with a minimum plasma 
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temperature of k B T ~ 50 keV and with the variations in the op- 
tical and UV light curves being due to thermal reprocessing. 
We could have chosen smaller values of the accretion rates in 
order to satisfy the limits imposed by M cr it, but that would ne- 
cessitate increasing the plasma temperatures in order to have 
larger thermal Comptonization luminosities, possibly violating 
our implicit constraint of kT < m e c 2 for NGC 3516 and wors- 
ening the situation for NGC 7469. Furthermore, the lower ac- 
cretion rates would require much greater radiative efficiencies, 
and for NGC 3516, the efficiency is already a fair fraction of 
unity. We could trade off reductions in the intrinsic disk lumi- 
nosity with a larger anisotropy factor, but values of £ > 1 would 
likely be required. 

However, for NGC 7469, our analysis really has only con- 
strained the quantity MM ^ 1.5 x 10 6 M 0 2 yr' 1 . If the cen- 
tral black hole mass were M ~ 10 8 M 0 , then we would have 
M ~ 0.03Af Edd which would just satisfy the accretion rate limit 
for advection for a ~ 0.2, assuming a Kerr metric. Since the re- 
verberation mapping estimates of the black hole mass in NGC 
7469 are given by M oc RblrVblr/^, where /?blr and v B lr are 
(model-dependent) length and velocity scales which character- 
ize the broad line region, these quantities would have to be each 
underestimated by factors of ~ 10 1 / 3 ~ 2 in order to account 
for an order-of-magnitude discrepancy in the central black hole 
mass (see Peterson & Wandel 2000). For a Schwarzschild met- 
ric, a central black hole mass of M ~ 3 x 10 7 M 0 would yield 
an accretion rate M ~ 0.04Af Edd . In this case, Rblr and v BE r 
would need only each be underestimated by a factor ~ 1.4. For 
M = 3 x 10 7 Mq, we have r^n/^g ~ 60 which would still be 
consistent with the inferred appearance of a narrow Fe Ka line 
in the RXTE/PCA spectra. 

By contrast, the required accretion rate for NGC 3516, as- 
suming M = 2 x 1O 7 M 0 and a Kerr metric, is M ~ 0.09 Af^. 
This too nominally exceeds M cr i t ; and since we have assumed 
r min = 6r gJ there is not much freedom here to increase the central 
black hole mass in order to increase M Edd and decrease M . If we 
assume r^ - 2 r gy which is still consistent with the Fe Ka emis- 
sion line fits of Nandra et al. (1999), the black hole mass could 
be larger by a factor ~ 3, and we would have M ~ 0.03A/ Edd . In 
any case, since the inner hot accretion flow we infer from our 
model is more “disk-like” than that of NGC 7469, a larger ac- 
cretion rate relative to Af Edd for NGC 3516 may nonetheless be 
consistent in some sense with the “strong ADAF proposal” of 
Narayan & Yi ( 1995; see also Narayan, Mahadevan, & Quataert 
1998), which states that NGC 3516 should then be less likely to 
have an ADAF-like flow than NGC 7469. Clearly, these machi- 
nations illustrate the speculative nature and uncertainties in- 
volved in extracting physical significance from our simple spec- 
tral fits based on the radiative processes. Regardless, the more 
plausible of the above suggestions are not wildly out of line 
with the observations or theoretical expectations. This perhaps 
indicates that this sort of modeling is on the cusp of actually 
providing useful constraints on the dynamical models. 

Numerous physical considerations which have not been in- 
cluded in these calculations may affect our results. In the case 
of a flattened Comptonizing region, such as we found for NGC 
3516, the mean Thomson depth through the plasma will vary 
significantly depending on the scattering order. Along with the 
non-uniform illumination of the Comptonizing region by the 
disk photons, this would lead us to expect to see features in 
the X-ray spectrum similar to the anisotropy breaks found by 
Stern et al. (1995) in their thermal Comptonization models of 


disk coronae. If the Comptonizing regions are threaded by suf- 
ficiently strong magnetic fields, then thermal synchrotron radi- 
ation can make a large contribution to the seed photon flux and 
thus reduce the importance of the disk radiation in determining 
the energy balance of the plasma as in the ADAF model of Esin 
et al. (2001). However, Wardzinski & Zdziarski (2000) have 
argued that for Seyfert Is the disk thermal emission should be 
a much more important source of seed photons than thermal 
synchrotron emission except for very low luminosity sources 
which have L/L Edd < 10 -4 . Finally, as we have already dis- 
cussed, enhanced illumination of the disk due to photons fol- 
lowing curved geodesics or anisotropy in the scattered radia- 
tion will clearly alter both the thermally reprocessed disk flux 
as well as the Compton reflection strength. 

Despite these deficiencies, the relative success of this model 
in reproducing the optical, UV and X-ray flux and spectral vari- 
ability observed from NGC 3516 and NGC 7469 underscores 
the crucial point that any model which attempts to describe the 
emission in these objects must at least take into account the 
roles that feedback and energy balance play in determining the 
emission across a wide range of energies. As corollary to this, 
observations of these objects would do well to consist of simul- 
taneous broad band spectral coverage, from the optical to the 
soft gamma-rays, in order to maximize the constraints on this 
and future, more sophisticated radiative and dynamical models. 
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APPENDIX 

FORMULAE FOR COMPUTING DISK INCIDENT FLUXES AND SEED PHOTON LUMINOSITIES 

Deriving these formulae from the specified geometry is straight-forward. We present the following expressions for completeness. 
Here all lengths are in units of r m iny i.e., r-rf /"min, etc.. 
r s > 1 (Spherical case, same as ZLS): 

C*max = 7r/2, f < f s , (Al) 

= sin~ 1 (r,/r), r>r s , (A2) 


^max ~ i f ^ 


= COS 


-1 


(P 2_p2) 1/2 


rcosa 


r>r 5 , 


l = cos a(r cos <f>+[r]( 1 + tan 2 a) - r 2 (sin 2 <f> + tan 2 a)] l/2 ), r < r s , 
= 2cosa[r?(l +tan 2 a)-r 2 (sin 2 0 + tan 2 a)] l A r > r s . 


r s < 1 (Ellipsoidal case): 


£*max “ tatl 


0max — COS 


L(^-i) 1/2 J 

[(r 2 - lXr^ + ta^a)] 1 / 2 


/ = 


2r 3 [(rj + tan 2 a) - /^(tan 2 a + rj sin 2 <j>)) ‘/ 2 
cosa(rf + tan 2 a) 


(A3) 

(A4) 

(A5) 

(A6) 

(A7) 

(A8) 

(A9) 
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Fig. A 1 Model geometries of the Comptonizing plasma. For r s > r mn , the plasma region is a sphere, and for r s < n it is an ellipsoid with semi-major axis 
r niin and semi-minor axis r s . A sample ray is shown with a path length / through the ellipsoidal Comptonizing region and making an angle a with respect to the disk 


pla 
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Fig. A3. — Model fluxes and fitted parameters for NGC 3516. The upper left plot shows the 2-10 keV X-ray light curve (solid curve). The dotted vertical lines 
are the boundaries of the epochs used to determine the X-ray spectral indices, and the horizontal line segments are the mean fluxes that were used to compute L x 
(see Eqs. 12 and 13). The middle left plot shows the optical continuum light curve (solid curve) averaged over the fluxes at 3590A, 4235A, and 5510A and scaled 
to unity. Also shown are the mean values for each epoch (line segments), and the fitted values for Mi = 1 (diamonds), 2 (triangles), and 3 (squares). The lower right 
plot shows the thermal Comptonization luminosities compared to the available accretion power for Mi =2 (dashed line) and Mi = 3 (dot-dashed line). The value for 
Mi - 1 is Tuntr — Tjisk = 3.36 x 10 44 erg s' 1 (see Table A2). 
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FlG. A4. — Model spectral energy distributions for NGC 3516 compared to the observed spectra. The SEDs are show for M 1 - 1 (solid curve), 2 (dashed), 3 
(dot-dashed). The diamonds are the de-reddened optical fluxes from Edelson et al. (2000) which have been used in the fits to the disk thermal spectra. The triangle 
in each plot is the de-reddened non-simultaneous flux at 1360A. The thick line segments are the 2-10 keV power-law continua which we have fit to the archival 
RXTE data. 
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Fig. a 5. — Model light curves and fitted parameters for NGC 7469 using anisotropy parameter values £ = 1 (filled diamonds), corresponding to isotropic thermal 
Compton emission, and £ = 1.5 (plus-signs); see Eq. 17. In the middle left plot, the measured fluxes are the open diamonds joined by line segments. The upper 
set of curves are the 13 15 A fluxes, and the lower set of curves are the 4865A fluxes. The dashed line in the lower left panel is the disk inner truncation radius, 
'•min = 3 x 10 14 cm. 
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Fig. A6. — Model light curves and fitted parameters for NGC 7469 using £ — 1.5 and allowing the X-ray spectral indices for each epoch to vary within the 
measured \-<7 bounds. The open diamonds are the measured or best-fit values determined from the data while the filled diamonds are the model values or parameters 
found from our spectral fits. 
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Fig. A7. A subset of the model spectral energy distributions for NGC 7469. The solid and dashed curves correspond to £ = 1 and £ = 1.5, respectively, using 

fixed values for the X-ray spectral indices, T. The dot-dashed curves are the SEDs found by allowing the T-values to vary within the 1-cr bounds and using £ = 1.5. 
The open diamonds are the starlight subtracted, de-reddened optical and UV fluxes from 1996 monitoring campaign The thick line segments are the 2-10 keV 
power-law continua from the fit parameters given in Nandra et al (2000). 
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p!G. A9. Compton reflection fraction versus X-ray spectral index for NGC 3516 (upper figure) and NGC 7469 (lower). The filled symbols are for R l= o = 

i/2Lx.zpp which is the implied reflection fraction for a face-on geometry, while the open symbols are for /?zls = Ld/(Lx ~ Id) which is the mean reflection fraction 
defined by ZLS. The reflection fractions and spectral indices obtained from fitting the pexrav model to the RXTE/PCA data are shown with the error bars. 





Table A 1 

Optical and X-ray data for the 1998 April Observations of NGC 3516 


start time 
(JD - 2450000.5) 

duration 

(ks) 

^3590 

(10- M 

^4235 

l erg cm'V 

^5510 

■A-') 

^2-10 

(i0~ n ergcnT 2 s -1 ) 

r 

916.692 

30.7 

1.8052 

0.9536 

0.6294 

6.75 

1.693 ±0.024 

917.048 

44.2 

1.8016 

0.9518 

0.6282 

5.53 

1.636±0.015 

917.559 

46.2 

1.8182 

0.9606 

0.6340 

6.29 

1.671 ±0.018 

918.094 

25.9 

1.8271 

0.9652 

0.6371 

5.81 

1.660 ±0.020 

918.393 

31.6 

1.8278 

0.9656 

0.6373 

6.03 

1.657 ±0.017 

918.759 

74.6 

1.7945 

0.9480 

0.6257 

4.94 

1.628 ±0.015 


Table A2 

Fixed Model Parameters, Luminosities, and Radiative Efficiencies for NGC 3516 


parameter 


value 


cos/ 


0.820 


dt (10 26 cm) 


1.22 


M (10 7 M o ) 

1 

2 

3 

Fmin (10 I3 cm) 

0.89 

1.78 

2.67 

M (10- 2 M Q yr‘) 

1.86 

0.946 

0.633 

^(lO^ergs- 1 ) 

5.26 

2.67 

1.79 

irfUk (10 44 erg s' 1 ) 

1.90 

0.97 

0.65 

Ai.avg (lO^erg s" 1 ) 

1.05 

0.75 

0.72 

( At,avg 4* Ltjjsg)/ Mc“ 

0.28 

0.32 

0.38 

( A*,avg 3" Ldisk)/^Edd 

0.23 

0.07 

0.03 


Table A3 

Fixed Model Parameters, Luminosities, and Radiative Efficiencies for NGC 7469 


parameter 


value 


cos i 


0.866 


di (10 26 cm) 


2.06 


M (10 7 M 0 ) 


1 


r min (10 13 Cm) 


*30. 


£ (anisotropy) 

1 

1.5 

1.5 

T (allowed variation) 

fixed 

± 1 cr 

M (Mgyr -1 ) 

0.159 

0.116 

0.150 

Una (10 44 erg s’ 1 ) 

45.0 

32.8 

42.5 

Ldi sk (lO^ergs -1 ) 

0.63 

0.46 

0.60 

At,av g (lO^erg s’ 1 ) 

1.43 

1.93 

1.47 

(^.r,avg 3" Ldjsk)/ 

0.023 

0.036 

0.024 

(Ai.avg 3" i'disk)/^Edd 

0.16 

0.19 

0.16 


